function J= computerEssentialJacobian(x,F)

f11=F(1,1);
f12=F(1,2);
f13=F(1,3);
f21=F(2,1);
f22=F(2,2);
f23=F(2,3);
f31=F(3,1);
f32=F(3,2);
f33=F(3,3);

fc=x(1,1);
ux=x(1,2);
vy=x(1,3);



J =[
 
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (8*fc^3*f11^2+8*fc^3*f12^2+(2*ux*f11+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(ux*f11+f12*vy+f13))*fc^2*f11+2*(2*fc^4*f11^2+2*fc^4*f12^2+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(8*fc^3*f11*f21+8*fc^3*f12*f22+(2*ux*f11+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+vy*f22+f23))*fc^2*f21+2*(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+(2*ux*f11+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f11+vy*f21+f31)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*fc^2*f11-2*(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f11,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f11)*fc^2*f11+(2*fc*f11*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f21)*fc^2*f21+(2*fc^3*f11^2+2*fc^3*f12^2+2*fc*f11*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f11*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*fc^2*f11,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      (2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f12)*fc^2*f11+(2*fc*f12*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f22)*fc^2*f21+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f12*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f21*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*fc^2*f11]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (8*fc^3*f11^2+8*fc^3*f12^2+(2*ux*f11+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(ux*f11+f12*vy+f13))*fc^2*f12+2*(2*fc^4*f11^2+2*fc^4*f12^2+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(8*fc^3*f11*f21+8*fc^3*f12*f22+(2*ux*f11+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+vy*f22+f23))*fc^2*f22+2*(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+(2*ux*f11+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f12+vy*f22+f32)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*fc^2*f12-2*(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f12,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f11)*fc^2*f12+(2*fc*f11*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f21)*fc^2*f22+(2*fc^3*f11^2+2*fc^3*f12^2+2*fc*f11*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f12*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*fc^2*f12,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      (2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f12)*fc^2*f12+(2*fc*f12*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f22)*fc^2*f22+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f12*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f22*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*fc^2*f12]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           (8*fc^3*f11^2+8*fc^3*f12^2+(2*ux*f11+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(ux*f11+f12*vy+f13))*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11^2+2*fc^4*f12^2+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f11*ux+fc*f12*vy+fc*f13))*(ux*f11+f12*vy+f13)+(8*fc^3*f11*f21+8*fc^3*f12*f22+(2*ux*f11+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+vy*f22+f23))*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f21*ux+fc*f22*vy+fc*f23))*(f21*ux+vy*f22+f23)+(6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+(2*ux*f11+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*(fc*f11*ux+fc*f12*vy+fc*f13)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(ux*f11+f12*vy+f13),                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  (2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f11)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11^2+2*fc^4*f12^2+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(2*fc*f11*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f21)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(2*fc^3*f11^2+2*fc^3*f12^2+2*fc*f11*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(2*ux*f11+vy*f21+f31+f12*vy+f13)-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(fc*f11*ux+fc*f12*vy+fc*f13)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f11,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              (2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f12)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11^2+2*fc^4*f12^2+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(2*fc*f12*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*fc*f22)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f12*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+(2*fc*f11*ux+2*fc*f12*vy+2*fc*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(f21*ux+2*vy*f22+ux*f12+f32+f23)-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(fc*f11*ux+fc*f12*vy+fc*f13)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f12]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (8*fc^3*f11*f21+8*fc^3*f12*f22+(2*f21*ux+2*vy*f22+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(ux*f11+f12*vy+f13))*fc^2*f11+2*(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(8*fc^3*f21^2+8*fc^3*f22^2+(2*f21*ux+2*vy*f22+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+vy*f22+f23))*fc^2*f21+2*(2*fc^4*f21^2+2*fc^4*f22^2+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+2*vy*f22+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f11+vy*f21+f31)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*fc^2*f21-2*(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f21,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      (2*fc*f21*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f11)*fc^2*f11+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f21)*fc^2*f21+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f21*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f11*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*fc^2*f21,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (2*fc*f22*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f12)*fc^2*f11+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f22)*fc^2*f21+(2*fc^3*f21^2+2*fc^3*f22^2+2*fc*f22*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f11+vy*f21+f31)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f21*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*fc^2*f21]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (8*fc^3*f11*f21+8*fc^3*f12*f22+(2*f21*ux+2*vy*f22+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(ux*f11+f12*vy+f13))*fc^2*f12+2*(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(8*fc^3*f21^2+8*fc^3*f22^2+(2*f21*ux+2*vy*f22+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+vy*f22+f23))*fc^2*f22+2*(2*fc^4*f21^2+2*fc^4*f22^2+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+2*vy*f22+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f12+vy*f22+f32)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*fc^2*f22-2*(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f22,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      (2*fc*f21*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f11)*fc^2*f12+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f21)*fc^2*f22+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f21*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f12*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*fc^2*f22,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (2*fc*f22*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f12)*fc^2*f12+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f22)*fc^2*f22+(2*fc^3*f21^2+2*fc^3*f22^2+2*fc*f22*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f12+vy*f22+f32)*fc+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f22*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*fc^2*f22]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           (8*fc^3*f11*f21+8*fc^3*f12*f22+(2*f21*ux+2*vy*f22+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(ux*f11+f12*vy+f13))*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f11*ux+fc*f12*vy+fc*f13))*(ux*f11+f12*vy+f13)+(8*fc^3*f21^2+8*fc^3*f22^2+(2*f21*ux+2*vy*f22+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+vy*f22+f23))*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f21^2+2*fc^4*f22^2+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f21*ux+fc*f22*vy+fc*f23))*(f21*ux+vy*f22+f23)+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+2*vy*f22+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*(fc*f21*ux+fc*f22*vy+fc*f23)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(f21*ux+vy*f22+f23),                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              (2*fc*f21*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f11)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f21)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f21^2+2*fc^4*f22^2+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(2*fc^3*f11*f21+2*fc^3*f12*f22+2*fc*f21*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(2*ux*f11+vy*f21+f31+f12*vy+f13)-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(fc*f21*ux+fc*f22*vy+fc*f23)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f21,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  (2*fc*f22*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f12)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^4*f11*f21+2*fc^4*f12*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*fc*f22)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^4*f21^2+2*fc^4*f22^2+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(2*fc^3*f21^2+2*fc^3*f22^2+2*fc*f22*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+(2*fc*f21*ux+2*fc*f22*vy+2*fc*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(f21*ux+2*vy*f22+ux*f12+f32+f23)-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(fc*f21*ux+fc*f22*vy+fc*f23)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f22]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(ux*f11+f12*vy+f13))*fc^2*f11+2*(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+vy*f22+f23))*fc^2*f21+2*(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(4*(ux*f11+vy*f21+f31)^2*fc+4*(ux*f12+vy*f22+f32)^2*fc)*(ux*f11+vy*f21+f31)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f11+vy*f21+f31)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*(ux*f11+vy*f21+f31)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(ux*f11+vy*f21+f31),                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   (2*fc^3*f11^2+2*fc^3*f12^2+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f11)*fc^2*f11+(2*fc^3*f11*f21+2*fc^3*f12*f22+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f21)*fc^2*f21+(4*fc^2*f11*(ux*f11+vy*f21+f31)+4*fc^2*f12*(ux*f12+vy*f22+f32)+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f11+vy*f21+f31)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f11*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f11+vy*f21+f31)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f11,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   (2*fc^3*f11*f21+2*fc^3*f12*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f12)*fc^2*f11+(2*fc^3*f21^2+2*fc^3*f22^2+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f22)*fc^2*f21+(4*(ux*f11+vy*f21+f31)*fc^2*f21+4*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f11+vy*f21+f31)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f21*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f11+vy*f21+f31)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f21]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                          (6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(ux*f11+f12*vy+f13))*fc^2*f12+2*(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+vy*f22+f23))*fc^2*f22+2*(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(4*(ux*f11+vy*f21+f31)^2*fc+4*(ux*f12+vy*f22+f32)^2*fc)*(ux*f12+vy*f22+f32)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(ux*f12+vy*f22+f32)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*(ux*f12+vy*f22+f32)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(ux*f12+vy*f22+f32),                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   (2*fc^3*f11^2+2*fc^3*f12^2+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f11)*fc^2*f12+(2*fc^3*f11*f21+2*fc^3*f12*f22+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f21)*fc^2*f22+(4*fc^2*f11*(ux*f11+vy*f21+f31)+4*fc^2*f12*(ux*f12+vy*f22+f32)+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f12+vy*f22+f32)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f12*fc-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*(ux*f12+vy*f22+f32)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f12,                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   (2*fc^3*f11*f21+2*fc^3*f12*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f12)*fc^2*f12+(2*fc^3*f21^2+2*fc^3*f22^2+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f22)*fc^2*f22+(4*(ux*f11+vy*f21+f31)*fc^2*f21+4*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f12+vy*f22+f32)*fc+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*f22*fc-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*(ux*f12+vy*f22+f32)*fc-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*fc*f22]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      (6*fc^2*f11*(ux*f11+vy*f21+f31)+6*fc^2*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(ux*f11+f12*vy+f13))*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f11*ux+fc*f12*vy+fc*f13))*(ux*f11+f12*vy+f13)+(6*(ux*f11+vy*f21+f31)*fc^2*f21+6*(ux*f12+vy*f22+f32)*fc^2*f22+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+vy*f22+f23))*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f21*ux+fc*f22*vy+fc*f23))*(f21*ux+vy*f22+f23)+(4*(ux*f11+vy*f21+f31)^2*fc+4*(ux*f12+vy*f22+f32)^2*fc)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)-(4*fc^3*f11^2+4*fc^3*f12^2+2*(fc*f11*ux+fc*f12*vy+fc*f13)*(ux*f11+f12*vy+f13)+4*fc^3*f21^2+4*fc^3*f22^2+2*(fc*f21*ux+fc*f22*vy+fc*f23)*(f21*ux+vy*f22+f23)+2*(ux*f11+vy*f21+f31)^2*fc+2*(ux*f12+vy*f22+f32)^2*fc)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33), (2*fc^3*f11^2+2*fc^3*f12^2+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f11)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f11+(2*fc^3*f11*f21+2*fc^3*f12*f22+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f21)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f21+(4*fc^2*f11*(ux*f11+vy*f21+f31)+4*fc^2*f12*(ux*f12+vy*f22+f32)+(4*ux*f11+2*vy*f21+2*f31+2*f12*vy+2*f13)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(2*ux*f11+vy*f21+f31+f12*vy+f13)-(2*fc*f11*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f21+2*fc^2*f11*(ux*f11+vy*f21+f31)+2*fc^2*f12*(ux*f12+vy*f22+f32)+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(2*ux*f11+vy*f21+f31+f12*vy+f13))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(2*ux*f11+vy*f21+f31+f12*vy+f13), (2*fc^3*f11*f21+2*fc^3*f12*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f11*ux+fc*f12*vy+fc*f13)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f12)*(fc*f11*ux+fc*f12*vy+fc*f13)+(2*fc^3*f11*(ux*f11+vy*f21+f31)+2*fc^3*f12*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f11*ux+fc*f12*vy+fc*f13))*fc*f12+(2*fc^3*f21^2+2*fc^3*f22^2+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*(fc*f21*ux+fc*f22*vy+fc*f23)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*fc*f22)*(fc*f21*ux+fc*f22*vy+fc*f23)+(2*fc^3*f21*(ux*f11+vy*f21+f31)+2*fc^3*f22*(ux*f12+vy*f22+f32)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(fc*f21*ux+fc*f22*vy+fc*f23))*fc*f22+(4*(ux*f11+vy*f21+f31)*fc^2*f21+4*(ux*f12+vy*f22+f32)*fc^2*f22+(2*f21*ux+4*vy*f22+2*ux*f12+2*f32+2*f23)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)+(2*(ux*f11+vy*f21+f31)^2*fc^2+2*(ux*f12+vy*f22+f32)^2*fc^2+((2*ux*f11+2*vy*f21+2*f31)*ux+(2*ux*f12+2*vy*f22+2*f32)*vy+2*ux*f13+2*vy*f23+2*f33)*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33))*(f21*ux+2*vy*f22+ux*f12+f32+f23)-(2*fc*f12*(fc*f11*ux+fc*f12*vy+fc*f13)+2*(fc*f21*ux+fc*f22*vy+fc*f23)*fc*f22+2*(ux*f11+vy*f21+f31)*fc^2*f21+2*(ux*f12+vy*f22+f32)*fc^2*f22+2*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)*(f21*ux+2*vy*f22+ux*f12+f32+f23))*((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)-(fc^4*f11^2+fc^4*f12^2+(fc*f11*ux+fc*f12*vy+fc*f13)^2+fc^4*f21^2+fc^4*f22^2+(fc*f21*ux+fc*f22*vy+fc*f23)^2+(ux*f11+vy*f21+f31)^2*fc^2+(ux*f12+vy*f22+f32)^2*fc^2+((ux*f11+vy*f21+f31)*ux+(ux*f12+vy*f22+f32)*vy+ux*f13+vy*f23+f33)^2)*(f21*ux+2*vy*f22+ux*f12+f32+f23)]
%[ 4*fc^3*f11*f22*f33-4*fc^3*f11*f23*f32-4*fc^3*f21*f12*f33+4*fc^3*f21*f13*f32+4*f31*fc^3*f12*f23-4*f31*fc^3*f22*f13,                                                                                                                 0,                                                                                                                 0]
 
 
];









end
